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Recently-performed research on developing simulation methods has /37 * 
made it possible to construct a general model of the biosphere Involving 
socio-economic, physico-chemical, and ecological processes. In following 
the basic concepts in models of the D. Forrester [1971] type, we shall 
consider the cause-and-effect links existing in the biosphere between 
such elements as the energy of solar radiation (E), the concentration 
of CO 2 (C) and O 2 (0), atmospheric turbidity (B), temperature (T), 
population (g, G), humus (S^, S^), and the nekton (r) and phytoplankton 
(<$») in the world's oceans (fig. 1). We shall divide dry-land vegetation 
into three types, differing in productivity and participating with 
varying degrees of intensity in other ecological processes: forests 
(P L » » agricultural vegetation (P^, Q^), and other vegetation (P^, Qj-}. 

In order to calculate socio-economic heterogeneities existing in 
the biosphere, dry-land area (S = 0.73805 X 10^ km 2 ) is divided into 
two regions, in each of which the processes under study may proceed 
at differing rates. 

It is given that a solar radiation energy of E q = 1.94 cal/cm 2 min. 
enters the biosphere and is used by biospheric photosynthetic elements. 
Here the value of E (t ) reaching the Earth's surface is determined by /39 
the atmosphere’s turbidity and may be calculated using the formula 
E (t) = E (t ) *exp (-aB - 8), where the coefficient of atmospheric 
absorption of solar energy due to dust and clouds a = 11.643 X 10 , 

and the transparency index of the pure atmosphere 6 = 6.487. The 
dustiness of the atmosphere is determined by the number of atmospheric 
dust particles resulting from dust storms (5.4 X 10^ tons/year), volcanic 
eruptions, solid and liquid fuel combustion (4.3 X 10^ tons/year), and 
expulsion from various kinds of metallurgical and chemical works, cement 
plants, and other sources (15 X 10^ tons/year). 


^Numbers in the margin indicate pagination in the foreign text. 
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Fig. 1. Diagram of the cause- /38 
and-effect ties in the global 
biogeocenotic biosphere simula- 
tion. Designations are explained 
in the text. 

It is considered that in /39 
ocean and dry-land areas the 
coefficients of solar energy 
use are = 0.07, K p = = 0.5. 

The productivity of the ocean 
and land ecosystems is determined 
by the intensity of photosyn- 
thesis, which depends upon the 
discrepancy between optimal and 
actual illumination, the amount 
of fertilizers applied, and 
pollution in the corresponding 
environments. The quantity of 
applicable fertilizers is deter- 
mined by the proportion of expend- 
able mineral resources (U , U_). Phytomass production increases in propor- 
tlon to the Increase In atmospheric CC^ concentration until C* = 0.2$ and 
decreases in proportion to the increase in 0^ content to approximately 
0* = 21$, so that 

(L = P,Q) 

As in the work of M. E. Vinogradov et al. (1973), we shall write 
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differential equations describing the change in phytomass (In tons/km ); 



where , Rp, and Rq are the rates of Increase in <f>, R, and Q, subject /40 
to illumination, the concentrations of C0 p and 0^, and environmental 
pollution. R d and R n also depend upon the distribution of sub-forest 







areas (S R ) and agricultural vegetation (S^). The formula for Rp, for 
example, is: 

Rp s KpPApEE"p‘exp[m 2 (<-EEp )][i-exp(-t 2 P)]> 

Where ; A -rS». |, ‘8 L l " *S V"B«' , .',(,.. si' 1 - Sx ! 'L',c. .■>] 

s, s, , ,3 l !■ 
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Here Ei is the optimal illumination for photosynthesis P (<ccal/m /day); 

/ -i \ W -i \ / -i \ 


B^ 1 ^, and b£‘ l; are P/B coefficients for L, X, and E, respectively; 


,(D 


M t 


, and Mq are the rates 


Ajjl is the coefficient of proportionality; M^, .-.p 
of <t> * P, and Q desiccation; t^, tp, and t^ are values for energy exchange 
expenditure; is the coefficient reflecting the value for the quota 
and proportion of the producer "j" in the nutritional allowance of 
consumer "i" (i,j = g, <p , P, Q, f, G, F). 


We shall assume that nekton is harvested by the regions at intensities 
of X g (t) and X G (t), maximum P/B is the coefficient r equal to K r , and the 
limitation in growth of r because of pollutions c and Z in regions 1 and 
II are described by the function ip = exp[-dr (pc+0Z ) ] , where p and 6 are 

proportions of all pollutions falling from regions 1 and II into the 
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ocean. Thus the change in biomass r (tons/km ) may, as ir the work of 
M. E. Vinogradov et al. (1973), be described by the equation: 

dt 1 ' (pr + Ag + /&) r ~ r r » 

where 

R* = K r r[l-exp( _ K r £p 1 J ) )] 1 ^ > 

y is the instantaneous mortality rate, and t is the value for energy 
exchange expenditure. 


We shall describe the population growth in both regions by calcu- 
lating the relationship of birth and mortality rates to the nutrition 
equation F Ra = V /a, the material standard of living /J}1 

'• ~ r i > - Q . rr _ rr \r f: j c W 0 l ' 1 
" j,; 'il '-'Ri Srit/E R u J , 

environmental pollution, atmospheric C0 2 and 0 o content, and population 
density. Here is the principal (basic resources) of the i region; 

S , U D . , and U„ . are the proportions of principal investments for the 
i region in the development of agriculture, renewal of mineral resources, 
and environmental conservation, respectively. 
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The equation for G is similar in form. 

Change in animal food in both regions is determined by the rate of 
animal growth, which depends upon the gaseous makeup of the atmosphere, 
the availability of a vegetable diet for animals, mortality rate, rate 
of energy exchange with the environment, and consumption by the populace. 
The equation for f, for example, has the following form: 

dt = f U,i " K gi/ R g v g < * 

where /42 

V f = K fp p + K fr IV * and v the ainou nt of harvestable nekton needed by 
the animals. 


In order to calculate the physiological effect of 0 o and concen- 
trations in living organisms, we shall cite the following relationship 


of the expenditure and energy exchange of element "a" 


t", and 
a’ 


t aiC + tao , C>C a i 
too, C 6 [o, c a ] ; 


f t ao, 0>0a i 


f-" -0 ff" -t n c h ' 
t a o onv L fl[i ^aoJ 


These relationships reflect the increased respiratory expenditure of 


animals and humans, with a CO^ concentration elevated beyond the 
threshhold of C and an 0 o concentration below 0. 

B fL a. 
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We shall descibe the effect of human agricultural activity on 
the environment by calculating the amount of energy required for popu- 
lation respiration, mineral resource utilization, and the generation* 
of pollution. We shall calculate possible ways to manage these elements 
which will make it possible to prevent environmental degradation by 
organizing the intelligent use and conservation of resources, and by 
imposing closed production cycles which fully utilize wastes and pollu- 
tion. Without detailing methods for achieving control, we shall propose 
that the following time functions are characteristic for the scientific 
and technological progress of both regions: (K^) is the pollution 

generation per individual person (a characteristic of the standard of 
living and the technology of production in a society); T* (TjJ) is the 

pollution resorption rate index; (M„) is the intensity of non-renewable 

M \ E « 

resource expenditure; TJL/ is the time interval necessary for transfer 

(i (i) 

to new resources in the i region; t£ and tg are the times needed by 
the i region for maximum incorporation of all lands suitable for culti- 
vation and achievement of the maximum possible productivity for 
agricultural cultivation. /**3 

We may write the following equation to describe the process of 
pollution generation and utilization fist region): 

g.K^tJgZvg-^j-UsgV^Cgg , 

where Z Vg 'Z?; x [^exo(-C^V ft9 )] , C 59 =C; g + C; 3 exp[-n Sfr Ct-t 0 )], 

T$(t) a T£ [t)[az g + a.' 2g Z r 9 9 ] • 

We shall describe the change in biospheric gaseous makeup by calcu- 
lating the following natural and anthropogenic processes. Assume that 
Og and C0 2 exchange between the atmosphere and ocean is described by 
L. Makht's (1971) model, the source of 0 2 output on dry land is the 
phytomass, 0 2 consumption occurs during respiration of element "a” at 
a rate of y a * v t a Wa and takes place in the resource combustion 
process at rates of 2> gQ and 2> GQ per person in regions 1 and II, respec- 
tively. Carbon Is assimilated by the plants „.f region 1 (II), in the 
form of C0 2 from the atmosphere, at rates of 0£ c R p (®qc R Q^ by forests, 
0 XP R P (®xQ R Q^ by a sricultural vegetation, and 6£ c R p (®qc r q^ b y other 
types of plants; it is liberated in the process of element "a" respiration 
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at a rate of t a“a ( given off during reBource combuatlon at a 

rate of t g( , (»„„>, and by the decomposition of dead vegetation at a rate 
of v s . 

We shall describe atmospheric turbidltv with the following equation: 
f? * N.S ♦ W,2 * W, t 9 c 9 • N, 6»c 5 - & *f» & 
where and (N2 and N^) are the amounts of pollution and smoke 
generated by region 1 (II) into the atmosphere, pg i§: the rate Of cloud 
cover alteration due to temperature fluctuations, T Q is the rate of 
natural clarification of the atmosphere due to dust settling. We shall /M 
present the following relationship between temperature, illumination, and 
alterations In biosphere gaseous makeup through the use of models such 
as that described by M. I. Budyko (1971). 

The remaining model equations, reflecting the dynamics of mineral 
resource alteration M g (M Q ), the principal V, (V 2 ) \nd capital investment 
(S 2 ) in agriculture investment for region 1, have the following 
form (Gelovani et al., 1976): 

- m 9 (t) R„ 9 g + V<U t9 G , 

ji'' kvgV M9 9 ‘V<T V9 , 

Sag' Sg ]T Sg , 

where R Mg = a R9 en(<*M sg ), C* 9 = C*R 9 +C* 9 e Uc9 ^ tm \ 

V«'«, S r ,'cxp(-5 S9 F„ ), 

V 6 ^ F ^ 13 


The model described above has been expressed in FORTRAN in the form 
of a program for the UVK M-*J030. Calculation of all Its components for 
100 years requires 1 hour of machine time, which makes it possible to 
carry out assorted experiments on the model. 

The prognosis for biosphere condition without any management and 
maintenance of rates of natural resource utilization, pollution generation, 
forest reduction, etc., is shown in fig. 2. Here, beginning in the year 
2050 , population density will fluctuate while maintaining a general 
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Fig. 2. Prognosis for biospheric component 
dynamics if contemporary rates of natural 
resource utilization are maintained. Explan- 
ation of the designations may be found in 
the text. Values for variables are in rel- 
ative units. Their values in 1970 are shown. 



Fig. 3. Prognosis for biospheric component 
dynamics if rates of influence on the bio- 
sphere and distribution of capital investment 
are decreased. Explanation in the text. 
Designations are the same as in fig. 2. 
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tendency to increase. By 2200 it will reach 1 85 persons/km 2 . The tem- 
perature of the lower atmosphere will increase 0.8°C, and this will lead ; 

to a rapid growth of vegetation (217 tons/km 2 in 2200). When t > 2200, 
existing food production rates will start to limit human population 
growth. Mineral energy sources will become the limiting factor after 
2300. Consequently, even if existing development rates for contemporary /46 
human society are maintained in the biosphere, acceptable conditions for 
human existence will continue for at least the next 250 years. During 
this period, mankind must first solve the problem of establishing an 
equilibrium with the environment and discover new energy sources. 

Now let's have a look at some hypothetical situations which might 
arise in the future. Figure 3 shows the results of a simulation which 
presupposed that by the year 2000 both regions will halve pollution 
generation rates, non-renewable resource expenditure will be reduced 50*, 
nekton harvesting will drop 10 %, agricultural capital will increase up 
to 45*, and capital investment U Rg * U R(J « 10* and U ?g * U ZQ “5*. By 
2000, 80 * of the land area suitable for farming will be utilized, animal 
husbandry productivity and P/B — the coefficient of agricultural vege- 
tation — increase 1.5 and 4 times, respectively, relative to 1970 levels, 
and values for T«, T*, and T Q decrease 30*. 

In this case, apparently, the system enters a quasi-stationary mode 

2 

where population density fluctuates from 50-200 persons/km . Either C0 2 
or food become the limiting factor at various stages. Non-renerable 
resources become limited after the year 2300. If mankind succeeds In 
switching to a new level of resource utilization by then, the "catastrophe" 
will be averted. 

Model calculations indicate that the worst gaseous conditions in 
the atmosphere may set in by 2070, if coordinated inter-region management 
is not achieved. CO^ concentration would exceed 0.0748*. 
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Fig. 4. Changes in atmospheric COp 
concentration when there is a dis- 
crepancy in natural resource utili- 
ztion rates for two biospheric 
regions. Explanation in the text. 
Atmospheric COp concentrations are 
expressed in percents. 


The model we examined here, 
as preliminary calculations have 
demonstrated, is flexible enough 
to reflect the natural and anthro- 
pogenic processes in the biosphere. 
It may, therefore, be used to study 
assorted hypothetical situations 
in order to find adoptable actions 
for controlling these processes. /48 
Allowance for regionality in the 
model may make it possible to 
evaluate the role of the regions 
in the fate of the biosphere. 
Specifically, fig. 4 shows the 
dependence of atmospheric COp 
concentration on the relationship 


between initial regional component conditions, with the assumption that 
region 1 implements measures to conserve the environment and region II 
retains the same rate of affecting It. Only if 0 >_ 50 can region II 
independently adopt measures to maintain a COp concentration within 
0 . 04 % . 
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